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Introduction 


The  U.S.  Army  Corps  of  Engineers  (Corps)  is  responsible  for  maintaining 
and  constructing  coastal  stmctures  nationwide.  The  Corps  currently  maintains 
over  1,500  coastal  structures.  These  structures  include  rubble  structures  such  as 
breakwaters,  jetties,  and  revetments  as  well  as  sea  walls,  timber  cribs,  and  float¬ 
ing  breakwaters.  The  majority  of  these  stmctures  are  of  mbble-mound  constrac- 
tion.  Most  of  these  stmctures  provide  protection  for  navigation  projects.  They 
prevent  sediment  from  accumulating  in  inlets  and  provide  shelter  for  navigation 
and  for  harbors. 

The  basic  deterministic  methodology  used  for  designing  breakwaters  and 
similar  mbble-mound  coastal  stmctures  was  developed  several  decades  ago  and 
includes  a  combination  of  analytical  techniques,  empirical  formulae,  and  physi¬ 
cal  scale  models.  Figure  1  shows  the  primary  failure  mechanisms  that  are  com¬ 
monly  addressed  in  modem  mbble  stmcture  design.  The  majority  of  these  fail¬ 
ure  mechanisms  are  addressed  in  the  design  process  using  empirical  formulae. 
The  degree  of  wave  transmission  is  typically  the  primary  functional  criterion. 

Preliminary  design  includes  defining  project  requirements,  and  then  deter¬ 
mining  environmental  design  conditions,  such  as  design  storm  characteristics 
over  the  stmcture’ s  intended  life.  Stmcture  design  follows  with  specification  of 
stmcture  type  and  layout,  and  armor,  underlayer,  and  core  material.  The  stmc¬ 
ture  cross  section  is  determined  by  specifying  the  crest  height  and  width,  type  of 
toe,  and  seaward  and  landward  slopes.  The  crest  height  and  width  are  the  lowest 
values  that  will  still  provide  minimum  allowable  overtopping  rates.  The  type  of 
toe  is  determined  based  on  the  local  subgrade  material,  length  of  seaside  slope, 
and  type  of  stmcture  armor.  The  armor  is  sized  to  prevent  more  than  a  few  per¬ 
cent  of  the  units  from  being  displaced  in  the  event  of  a  design  level  storm.  The 
design  process  is  traditionally  deterministic,  although  some  elements  of  proba¬ 
bility  theory  have  been  utilized.  It  has  been  common  practice  to  specify  wave 
height  as  a  stochastic  parameter  associated  with  a  design  storm  with  a  given 
return  period.  But  many  more  stochastic  variables  must  be  included  in  order  to 
adequately  assess  the  reliability  of  a  coastal  stmcture  design. 

Reliability  methods  utilize  the  same  probability  theory  concepts  required  for 
specifying  a  design  wave  or  for  determining  confidence  limits  and  are  therefore 
quite  simple  to  apply.  They  are  readily  adaptable  to  a  wide  variety  of  coastal 
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Figure  1 .  Breakwater  failure  modes 

structure  design  and  evaluation  problems  and  provide  a  powerful  tool  for  ratio¬ 
nally  making  economic  compromises  that  are  always  necessary  in  civil  engineer¬ 
ing.  This  report  discusses  reliability  techniques  that  facilitate  inclusion  of  the 
random  nature  of  the  rubble-mound  response  in  the  design  process.  These  reli¬ 
ability  techniques  can  be  utilized  to  gain  insight  into  the  degree  of  uncertainty  in 
the  empirical  design  equations.  This  is  very  important  for  rubble  mounds  as 
their  response  varies  widely,  even  in  controlled  laboratory  conditions.  But  it 
must  be  kept  in  mind  that  these  methods  do  not  guarantee  a  more  reliable  struc¬ 
ture.  Limited  local  data,  the  lack  of  design  guidance  for  certain  failure  modes 
(e.g.  scour),  and  unknown  and  somewhat  unpredictable  construction  quality, 
among  other  things,  can  reduce  a  structure’s  reliability  far  below  the  design 
engineer’s  estimate.  Therefore,  while  reliability  methods  provide  a  powerful 
tool  for  gauging  and  comparing  the  uncertainty  of  certain  failure  modes  in  our 
designs,  the  actual  probability  of  failure  of  the  as-built  structure  may  be  signifi¬ 
cantly  different. 

Reliability  methods  can  be  a  powerful  tool  in  estimating  uncertainty  with 
respect  to  many  failure  modes.  Consider,  for  example,  armor  stability  design. 
The  preliminary  design  process  for  specifying  the  armor  layer  is  deterministic 
and  can  be  characterized  as  determining  a  design  load  and  then  defining  the 
appropriate  structural  capacity.  This  is  done  with  a  single  empirical  stability 
equation  which  incorporates  a  wave  height  statistic  from  an  assumed  distribution 
associated  with  a  particular  storm  return  period.  The  design  is  then  evaluated 
using  scaled  physical  models.  The  deterministic  design  method  has  several 
weaknesses,  as  follows: 

a.  Only  the  loading  variables  associated  with  the  wave  are  specified 
according  to  measured  distributions.  The  capacity  variables  are  typi- 
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cally  best-fit  mean  values  of  widely  scattered  quantities.  Confidence 
limits  are  rarely  specified  with  the  formulae  and  provide  little  help  to  the 
designer  for  determining  the  risk  or  reliability  in  a  given  design  or  in  the 
future  performance  of  an  existing  structure. 

b.  Laboratory  experiments  designed  to  validate  design  formulae  are  done 
under  highly  controlled  conditions  and  the  prototype  may  be  quite  dif¬ 
ferent  with  more  widely  varying  parameters. 

c.  The  equations  do  not  include  the  interaction  between  failure  modes,  such 
as  structural  response  and  hydraulic  stability  or  primary  armor  stability 
and  toe,  crown,  or  transitional  armor  stability.  They  typically  give  no 
indication  of  future  performance  of  a  damaged  structure. 

Given  the  uncertainties  in  deterministic  breakwater  design  using  empirical 
equations,  there  is  need  to  incorporate  objective  evaluation  of  the  performance 
risks  of  either  new  or  existing  structures.  Coastal  structures  can  be  allowed  to 
sustain  considerably  more  damage  than  land-based  structures  because  human 
lives  are  typically  not  in  jeopardy.  Engineering  reliability  methods  can  be  used 
to  compare  repair  scenarios  for  a  given  structure  or  to  rank  various  projects. 

With  repair  costs  included  in  the  economic  analysis,  the  most  economical  plan 
meeting  the  specified  reliability  can  be  selected.  This  evaluation  technique  pro¬ 
vides  more  economically  efficient  coastal  rubble  structure  designs  and  provides 
allowance  for  periodic  repairs. 

It  is  standard  practice  within  the  Corps  of  Engineers  to  use  reliability  meth¬ 
ods  to  compare  competing  engineering  alternatives.  But  for  coastal  structural 
engineering,  little  work  has  been  done  to  synthesize  data  collected  over  several 
decades  into  a  usable  format  for  reliability  analysis.  Carver  (1983)  performed  a 
rubble-mound  armor  stability  experiment  investigating  the  variability  in  stability 
due  to  the  natural  variations  in  stone  placement.  As  one  would  expect,  the  wave 
height  corresponding  to  the  no-damage  condition  varied  widely.  Using  these 
data,  it  is  shown  in  Chapter  3  of  this  report  that  the  variability  in  stability  can  be 
accounted  for  in  the  design  process. 

Typically  for  design,  all  dominant  failure  modes  must  be  identified  and  a 
fault  tree  constructed  for  each  design  alternative.  A  fault  tree  is  a  flowchart  of 
failure  modes.  An  example  of  a  fault  tree  for  a  breakwater  analysis  is  shown  in 
Figure  2.  The  references  herein  provide  guidance  on  establishing  a  fault  tree. 
Moritz  et  al.  (1994)  performed  a  reliability  analysis  of  engineering  alternatives 
for  a  major  rehabilitation  of  the  Bums  Waterway  Harbor,  Indiana,  breakwater. 
They  describe  the  fault  tree  and  evaluation  of  alternatives.  Although  the  reliabil¬ 
ity  calculations  are  not  shown,  they  provide  an  overview  of  the  entire  alterna¬ 
tives  analysis  utilizing  reliability  methods  and  an  economic  analysis. 

Although  several  reliability  analysis  methods  have  been  used  for  coastal 
structures,  no  rigorous  comparative  analysis  of  the  various  methods  has  been 
published.  Meadowcroft  et  al.  (1996)  compared  the  Level  III  Monte-Carlo  tech¬ 
nique  to  a  Level  II  mean-value  approach.  Although  details  of  the  study  were  not 
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described,  they  found  that  the  Level  II  approach  “provided  reasonable  agree¬ 
ment”  to  the  more  laborious  Level  III  approach.  The  mean  value  approach  is 
described  in  the  following  two  chapters  as  is  a  more  accurate  Level  II  method 
termed  the  “design  point  approach.” 

This  report  presents  an  introduction  to  and  comparison  of  various  Level  II 
reliability  methods  as  applied  to  coastal  rubble  structure  design.  Several  compu¬ 
tational  techniques  of  varying  degrees  of  accuracy  are  presented.  Examples  of 
computational  methods  are  shown  in  tabular  format  and  FORTRAN  programs. 


Figure  2.  Partial  fault  tree  for  breakwater  analysis 
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2  Engineering  Reliabiiity 
Methods 


Reliability  Methods  Overview 

Engineering  Technical  Letter  (ETL)  1 1 10-2-532  (U.S.  Army  Corps  of  Engi¬ 
neers  1992)  provides  a  review  of  reliability  methods  and  the  underlying  proba¬ 
bility  theoiy.  Other  documents  that  provide  good  reference  are  Harr  (1987), 
Benjamin  and  Cornell  (1970),  Ang  and  Tang  (1975),  and  Thoft-Christensen  and 
Baker  (1982). 

Reliability  or  conversely,  probability  of  failure,  is  computed  from  the  proba¬ 
bility  distribution  of  a  limit  state  equation.  Several  reliability  definitions  that  are 
useful  follow. 

a.  Reliability.  Probability  that  limit  state  equation  will  be  greater  than  limit 
state. 

b.  Probability  of  Failure.  Probability  that  limit  state  equation  will  be  less 
than  limit  state. 

c.  Limit  State  Equation  (otherwise  known  as  the  failure  function  or  perfor¬ 
mance  function).  Equation  describing  the  engineering  performance  of 
interest  expressed  as  either  the  difference  between  capacity  and  demand 
(safety  margin)  or  ratio  of  capacity  to  demand  (safety  factor).  Some 
authors  use  the  words  resistance  and  load  rather  than  capacity  and  de¬ 
mand. 

d.  Safety  Factor.  Ratio  of  capacity  to  demand. 

e.  Safety  Margin.  Difference  between  capacity  and  demand. 

f.  Limit  State.  Level  of  performance  for  which  capacity  equals  demand 
(safety  factor  =  1  and  safety  margin  =  0). 
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g.  Failure  Surface.  Surface  along  the  limit  state  described  by  the  limit 
state  equation. 

To  illustrate  these  definitions,  consider  breakwater  armor  stability  characterized 
by  the  Hudson  (1958)  equation. 


W  = 


liTo  (5^-1)^  cote 


(1) 


where 

W  =  weight  of  armor  unit 

Yr  =  specific  weight  of  armor  unit  material 

H  =  design  wave  height  at  the  structure  toe 

=  specific  gravity  of  armor  unit  material 

0  =  sea-side  angle  of  armor  slope  relative  to  the  horizontal 

K^)  =  tabulated  empirical  stability  coefficient 

For  this  equation,  /sT^,  is  defined  for  a  given  level  of  performance,  typically  the 
no-damage  condition  represented  by  less  than  2  percent,  by  count,  of  the  stones 
displaced  from  the  seaward  structure  face.  To  begin  the  reliability  analysis,  the 
equation  of  interest  is  usually  rearranged  to  form  a  safety  factor,  F=  C/D,  where 
C  represents  capacity  and  D  represents  demand.  In  this  case,  the  limit  state  is 
prescribed  when  the  safety  factor  is  equal  to  one.  Performance  is  satisfactory  if 
the  safety  factor  is  greater  than  one  and  unsatisfactory  when  the  safety  factor  is 
less  than  one.  For  a  safety  factor  approach,  the  Hudson  limit  state  equation 
could  be  formulated  as 


W^o(5^-l)3cot0 


(2) 


Another  way  to  express  the  limit  state  equation  is  as  a  safety  margin,  g  =  C  -  D. 
The  condition  under  which  this  margin  is  less  than  zero  prescribes  unsatisfactory 
performance,  the  condition  above  zero  satisfactory  performance,  and  the  condi¬ 
tion  at  zero  describes  the  limiting  state  of  performance.  The  Hudson  equation 
can  be  rewritten  in  the  form  of  a  safety  margin  as 


g  =  WK^iS^-lfcotd-y^H^ 


(3) 
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An  equivalent  safety  margin  can  also  be  formulated  as 


g  =  AD„  (Kj.cotey'^  -  H 


(4) 


where  D„  =  is  the  nominal  diameter  of  the  armor  unit  and  L=S^  - 1 . 

Because  some  or  all  of  the  variables  in  the  limit  state  equation  are 
nondeterministic,  the  limit  state  equation  is  also  nondeterministic.  The  process 
of  defining  the  reliability  requires  defining  probability  density  functions  (pdf)  of 
all  stochastic  variables  in  the  performance  function  and  combining  these  to  de¬ 
termine  the  overall  reliability  of  the  given  performance  function.  The  reliability 
is  computed,  using  probability  theory,  as  the  probability  that  the  performance 
function  will  exceed  the  limit  state.  For  the  Hudson  equation  characterized  by  a 
safety  factor  in  Equation  2,  the  reliability  is  the  probability  that  the  weight-domi¬ 
nated  capacity  parameter  in  the  numerator  will  exceed  the  wave  loading  parame¬ 
ter  in  the  denominator,  i.e.  the  probability  that  F  will  exceed  one. 

Several  reliability  estimation  methods  are  commonly  used  in  civil  engineer¬ 
ing  to  determine  the  reliability  for  a  particular  limit  state  equation  given  the  pdfs 
of  the  underlying  variables.  These  include  Level  I  methods,  where  a  coefficient 
is  used  in  the  design  equation  to  account  for  uncertainties  in  the  design.  This  is 
the  method  used  for  the  familiar  Load  and  Resistance  Factor  Design  (LRFD) 
structural  design  method.  Level  II  methods  provide  approximations  to  the  reli¬ 
ability,  assuming  the  limit  state  equation  to  be  normally  distributed  and  then 
converting  all  the  random  correlated  non-normally  distributed  variables  to  non- 
correlated  normally  distributed  variables,  or  assuming  a  mathematically  simpli¬ 
fied  form  of  the  failure  surface,  or  both.  Level  III  methods  utilize  the  actual 
distributions  of  the  random  variables  to  compute  the  reliability  of  the  limit  state 
equation.  Specific  techniques  for  the  Level  II  and  III  methods  are  discussed 
below. 

a.  Direct  Integration.  A  Level  HI  method  where  the  limit  state  reliability 
or  conversely,  the  probability  of  failure,  is  computed  by  integrating  or 
convolving  the  capacity  and  demand  pdfs. 

b.  Monte  Carlo  Simulations.  A  Level  III  method  where  the  limit  state 
reliability  is  approximated,  through  a  large  number  of  realizations  of 
capacity  and  demand  pdfs,  by  the  proportion  of  realizations  where  the 
limit  state  equation  is  greater  than  the  limit  state. 

c.  Taylor  Series  Approximations.  A  Level  II  method  where  the  limit  state 
equation  is  approximated  by  a  Taylor  Series  expansion  about  some  criti¬ 
cal  point  usually  truncated  to  first  or  second  order.  The  reliability  is 
computed  as  the  minimum  distance  between  the  failure  surface  and  zero. 

d.  N-Point  Estimates.  A  Level  II  method  where  the  pdf  of  the  limit  state 
equation  is  approximated  as  N  lumped  masses. 
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Taylor  series  and  N-Point  estimates  typically  provide  reasonable  accuracy 
without  having  to  perform  thousands  of  calculations.  They  are  particularly  use¬ 
ful  when  data  are  sparse  and  use  of  empirical  equations  derived  from  data  entails 
a  great  deal  of  uncertainty.  Such  is  typically  the  case  with  coastal  structure  de¬ 
sign.  Level  II  methods  also  provide  insight  into  the  sensitivity  of  reliability  to 
various  parameters  throughout  the  calculation  process  that  is  useful  for  the  engi¬ 
neer.  Several  popular  Level  II  methods  are  compared  herein.  More  accurate 
Level  III  methods  can  be  used  to  determine  the  reliability  if  the  density  functions 
of  all  nondeterministic  variables  are  known  or  assumed;  but  convolving  more 
than  two  density  functions  is  computationally  intensive  and  is  therefore  not  gen¬ 
erally  practical.  Also,  for  coastal  structure  engineering,  the  probability  distribu¬ 
tions  of  the  various  parameters  are  seldom  known  with  any  certainty. 


Reliability  Concepts 

As  stated  above,  the  first  step  to  estimate  the  reliability  of  a  design  is  to  es¬ 
tablish  a  limit  state  equation  in  one  of  two  forms: 


safety  margin:  g{x)  =  CiXff)  -  D{Xjf)  (5) 

safety  factor.  Fix)  =  CiXf.)/DiXp)  (6) 


such  that  g>0  or  F>1  represents  satisfactory  performance.  In  these  equations,  x 
represents  a  set  of  stochastic  variables  describing  geometry,  material  properties, 
and  loading  for  the  particular  limit  state,  C  is  a  function  of  the  set  of  stochastic 
capacity  variables  X(.,  emd  D  is  a  function  of  the  set  of  stochastic  demand  vari¬ 
ables  Xp.  Accordingly,  the  limit  state  function  is  also  stochastic.  Consequently  g 
and  F  can  be  described  by  pdf  s,  as  illustrated  in  Figure  3.  Note  that  gov  F  can 
be  a  function  of  many  variables.  The  failure  surface  is  defined  by  the  limit  state 
g=0  or  F=l.  If  the  limit  state  is  a  function  of  more  than  two  variables,  say  N 
variables,  an  imaginary  failure  "surface"  in  N-dimensional  hyperspace  can  be 
described. 

Reliability  with  respect  to  the  performance  function  g  can  be  computed  from 
the  relation  /?=P(g>0),  which  states  that  reliability  R  is  the  probability  of  satis¬ 
factory  performance  (g>0).  P  represents  the  exceedance  probability  distribution 
of  g.  Conversely,  the  probability  of  failure  with  respect  to  the  performance 
function  is  given  by  F^=  F(g<0).  Similarly,  for  a  safety  factor  representation 
^=^(^■>1)  and  F^=  F(F<1).  Figure  3  illustrates  several  of  these  points.  The 
region  under  the  pdf  to  the  left  of  the  origin  represents  the  probability  of  failure 
and  that  to  the  right  represents  the  reliability.  Note  that  the  tevm  failure  is  not 
used  to  describe  the  structure’s  performance  directly.  It  is  used  to  describe  the 
fact  that  the  performance  function  does  not  meet  the  design  criteria.  For  exam- 
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pie,  if  the  performance  function  for  armor  stability,  say  the  Hudson  equation, 
describes  the  design  condition  for  “2-percent  armor  displacement,”  then  the  term 
failure  would  describe  the  limiting  condition  of  2  percent  of  the  armor  being 
displaced.  This  contrasts  with  the  conceptual  picture  most  of  us  might  have  of 
breakwater  failure  where  the  rubble  mound  has  collapsed  or  been  breached. 


Often  the  limit  state  equation  pdf  is  not  completely  known  but  the  mean  and 
standard  deviation  are  known.  Based  on  the  principle  of  maximum  entropy,  the 
normal  distribution  can  be  assumed  if  only  the  first  two  moments  are  known 
(Harr  1987).  If  capacity  and  demand  functions  are  assumed  to  be  normally  dis¬ 
tributed,  then  the  limit  state  relation  will  also  be  normally  distributed.  If  capac¬ 
ity  and  demand  are  expressed  in  standard  normal  form 


C'  = 


C  -  Pc 


D' 


(7) 


then  the  limit  state  relation  (g  =  C  -  D  =  0)  becomes 

OcC  -  OjjD'  +  Pc  -  Po  =  0 


(8) 
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where  iJ^  and  are  the  mean  and  standard  deviation  of  the  capacity  function 
and  Pd  and  Oj,  are  the  mean  and  standard  deviation  of  the  demand  function.  Note 
g(C,D)  =  0  is  displaced  a  distance  p  from  the  origin  in  the  transformation  to 
g{C',D')  =  0.  This  distance  can  be  found  from  simple  geometry  (Hasofer  and 
Lind  1974)  as 


P 


(9) 


if  C  and  D  are  uncorrelated  (Figure  4).  P  is  termed  the  reliability  index.  Here  Pg 
and  Og  are  the  mean  and  standard  deviation  of  the  safety  margin  performance 
function.  If  C  and  D  are  correlated,  the  expression  becomes 


^PcD^C^D 


(10) 


where  is  the  correlation  coefficient  between  C  and  D  defined  as 


■'CD 


E[(C  -  \ic)(.L  -  p^)] 


(11) 


where  0^^  is  the  covariance  between  C  and  D  and  E[  ]  is  the  expected  value  of 
the  quantity  inside  the  brackets.  Note  that  P  is  simply  a  unit  standard  normal 
variate  along  z  =  (x  -  p,()/a,j  =  0.  The  reliability  can  therefore  be  approximated 
using  the  reliability  index  as 

R  =  P(g>0)  =  I  Pg(x)dx  =  a>(p)  (12 

0 


where  pg(x)  is  the  pdf  of  g(x)  and  <I>(P)  is  the  widely  tabulated  standard  normal 
distribution  function  evaluated  at  p.  The  definitions  for  F(x)  would  be  similar; 
but  the  log-normal  distribution  would  be  used  because  the  normal  distribution  is 
an  appropriate  approximation  for  sums  of  variables  while  the  log-normal  is  ap¬ 
propriate  for  products. 

The  reliability  index  is  graphically  defined  in  Figures  3  and  4.  It  can  be  seen 
that  the  index  represents  the  distance  from  the  expected  value  of  the  perfor¬ 
mance  function  to  the  failure  surface  in  units  of  standard  deviation.  Thus  P  =  2 
implies  that  the  mean  of  the  performance  function  lies  two  standard  deviations 
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above  the  limit  state.  The  reliability  index  is  then  a  measure  of  the  likelihood  of 
the  structure  to  perform  satisfactorily,  with  respect  to  the  performance  function 
in  question.  Reliability  increases  with  increasing  mean  and  with  decreasing 
standard  deviation  of  the  performance  function.  Table  1  shows  the  correspon¬ 
dence  between  the  reliability  and  the  reliability  index.  This  correspondence  will 
be  valuable  later  when  interpreting  applications  of  the  reliability  methods.  In 
Table  1,  4>  '  is  the  inverse  of  the  normal  distribution  function. 


Table  1 

Correspondence  Between  Reliability  Index  and  Reliability 

R=Pici>0)=P(F>1) 

3  =<D-’{fl) 

0.500 

0.0 

0.690 

0.5 

0.840 

1.0 

0.933 

1.5 

0.977 

2.0 

0.9987 

3.0 

0.999968 

4.0 

In  general,  g  will  be  a  nonlinear  function  of  the  random  variables  z.  The  case 
for  two  independent  random  variables  is  depicted  in  Figure  4.  Note  that  the 
transformation  into  a  normalized  coordinate  system  has  been  performed  where 
z  =  (^  -  |4)/a^.  The  reliability  index  is  the  minimum  distance  from  the  mean 
state  z  =  0  to  the  limit  state  g(z)  =  0.  The  closest  point  on  g{z)  =  0  is  known  as 
the  design  point  Zj. 


Figure  4.  Example  sketch  of  linear  and  nonlinear  limit  states 
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Often  in  computations,  the  limit  state  relation  is  approximated  as  a  linear 
failure  surface  using  a  Taylor  series  expansion  either  about  the  mean  point  or 
about  the  design  point.  An  important  characteristic  of  the  design  point  expan¬ 
sion  is  that  the  resulting  reliability  index  is  invariant  to  the  arbitrary  functional 
form  of  the  limit  state  expression;  whereas  a  mean  point  approximation  will  be 
dependent  on  the  functional  form  of  the  limit  state.  The  two  methods  will  be 
compared  in  the  next  section. 

If  the  random  variables  z  are  non-normal,  an  approximate  normal  distribution 
can  be  fitted  to  the  actual  distribution  such  that  the  density  and  distribution  are 
unchanged  at  the  design  point.  With  the  design  point  normal  variate  designated 
Zj  =  (Xj  -  |a,)/a,,  the  mean  and  standard  deviation  of  the  fitted  density  function 
are 


Px(Xd'>  (13) 


where  P  and  p  are  the  distribution  and  density  functions  of  random  variable  x,  $ 
and  (j)  are  standard  normal  distribution  and  density  functions,  and  $  ’  is  the  in¬ 
verse  of  the  standard  normal  distribution. 


Taylor  Series  Approximations  of  Reliability 

Taylor  series  mean  point  method 

Based  on  the  preceding  discussion  it  appears  that  a  simple  method  for  esti¬ 
mating  the  reliability  would  be  to  approximate  the  limit  state  either  at  the  mean 
or  design  point  on  the  failure  surface.  This  can  be  done  using  a  Taylor  series 
expansion  of  the  function  about  either  point.  The  mean  point  expansion  can  be 
done  directly;  but,  because  the  design  point  is  a  priori  unknown,  the  design  point 
solution  requires  an  iterative  approach.  Either  way,  g  is  first  linearized  about  the 
mean  point  x  =  0.  The  Taylor  series  approximation  of  the  safety  margin  about 
the  mean  point  is  given  by 


8  =  g(Kf)  +  Vg«(x  -  np 


(14) 


or  in  scalar  notation 


8  =  ^(l^f) 


(15) 
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where  =  is  used  to  denote  a  first  order  approximation,  Vg  •  is  the  dot  product 
operation  using  the  gradient  of  g,  partial  derivatives  are  evaluated  at  each  mean 
value,  and  N  is  the  number  of  random  variables  in  g.  The  general  expressions 
for  the  expected  value  and  variance  of  g  as  a  function  of  correlated  random  vari¬ 
ables  are 


=  Mg  =  +  Vg«(f  -  p-)]  =  g(p^-)  (16) 

V[g]  =ol^Vg  vm  '  (17) 


where  E  is  expected  value  and  V  is  variance  and  Vg^  is  the  transpose  of  Vg.  For 
N  variables,  the  variance  becomes 


N  AT 


El 

f=i  i=i 


^  a  ^ 
dx.dxj 


(18) 


and,  if  the  variables  are  uncorrelated,  the  variance  is 


dx. 


E 

j=i 


V 


(19) 


The  reliability  index  and  reliability  can  be  computed  from  the  mean  and  variance 
using  equations  10  and  12,  respectively. 

For  the  safety  factor  solution,  the  expected  value  and  variance  can  be  com¬ 
puted  using  the  log  of  the  mean  and  veuriance  of  the  safety  factor  relation  approx¬ 
imated  at  the  mean  or  design  points  according  to  Benjamin  and  Cornell  (1970) 
as 


Fmr  =  In  Ff  - 


(20) 


(f  ) 

2  ' 

+  1 

A 

/ 

(21) 


where  pp  =  F(p,j),  as  in  Equation  16,  and  Op  is  similar  to  Equation  18.  Then  the 
reliability  index  is  P  =  Pi„p/o,^  and  reliability  R  =  $(P). 
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Taylor  series  finite  difference  approximation 

The  Taylor  series  mean  point  approximation  is  typically  employed  using  a 
finite  difference  approximation  for  the  partial  derivatives  in  the  variance.  The 
calculation  is  often  performed  using  the  difference  between  the  limit  state  rela¬ 
tion  at  one  standard  deviation  above  the  mean  and  one  standard  deviation  below 
the  mean  as  follows 


dx.  2a. 


(22) 


This  equation  is  often  written  in  shorthand  as 

dg  ^  g"  -  g~ 

dx.  2aj 


(23) 


The  safety  factor  equations  would  be  similar.  The  mean  point  analysis  with 
finite  differencing  is  commonly  known  as  the  Taylor  Series  Finite  Difference 
(TSFD)  method.  The  strength  of  this  technique  of  reliability  estimation  is  that  it 
provides  a  closed  form  solution  that  can  be  done  easily  using  a  handheld  calcula¬ 
tor  or  a  spreadsheet  program.  The  weakness  of  the  mean  point  method  is  that  the 
solution  is  dependent  on  the  form  of  the  limit  state  equation,  i.e.  p  and  /?  can 
vary  depending  on  the  form  chosen.  Additionally  it  is  less  accurate  than  more 
rigorous  solutions.  Example  calculations  are  shown  in  Chapter  3. 

Taylor  series  design  point  method 

Mlakar  (1994)  proposed  a  more  accurate  iterative  solution  for  the  design 
point  definition  of  the  reliability  index.  Consider  Zj  =  XV  g(z)  as  an  iteratively 
improved  estimate  of  the  design  point,  where  A  is  an  arbitrary  constant.  Then 
substituting  into  Equation  14,  expressed  as  a  function  oft  yields  the  design 
point  on  the  linearized  failure  surface  which  is  closest  to  the  mean  as 


^  z-Vg(^  -  g(z) 
Vg(z)-Vg(z) 


yg(z) 


(24) 


Pd  =  Izrfl  now  represents  an  estimate  of  the  reliability  index.  This  solution  is  simi¬ 
lar  to  TSFD  but  more  generalized.  The  iterations  start  with  a  mean  point  ap¬ 
proximation.  At  each  iteration  step,  a  new  design  point  is  computed  from  Equa¬ 
tion  24  and  the  iterations  continue  until  convergence  is  achieved.  Typically  the 
partial  derivatives  are  numerically  approximated  using  a  finite  difference  formu¬ 
lation  such  as 
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dg  _  g(z+A)  -  gjz) 


(25) 


dz 


A 


where  A  is  some  small  increment  in  z.  Hereafter,  this  iterative  TSFD  method 
will  be  referred  to  as  ITSFD. 

A  general  FORTRAN  program  for  the  ITSFD  computation  is  shown  in  Ap¬ 
pendix  A.  The  main  program  comprises  the  reliability  engine,  while  the  specific 
limit  state  equation  is  computed  in  subroutine  LIMSTATE.  Subroutine 
DEP_MEAN  is  used  to  compute  a  range  of  one  variable  in  the  limit  state  equa¬ 
tion  so  a  range  of  reliabilities  can  be  computed  at  one  time.  An  example  break¬ 
water  limit  state  equation  is  encoded  in  Appendix  A  and  example  computations 
are  discussed  in  Chapter  3. 


Point  Estimate  Method 

Point  estimate  methods  (PEM),  as  described  in  Harr  (1987)  and  ETL  1 1 10-2- 
532  (USAGE  1992)  provide  invariant  approximations  to  the  reliability  index. 

The  N-point  estimates  proceed  in  a  similar  fashion  to  the  Taylor  series  approxi¬ 
mations  except,  rather  than  the  limit  state  relation  being  approximated,  the  prob¬ 
ability  distributions  are  approximated  at  N  points.  PEMs  have  the  advantage  that 
they  do  not  require  computation  of  derivatives  of  the  limit  state  equation  nor  do 
they  impose  assumptions  on  the  existence  and  continuity  of  the  first  few  deriva¬ 
tives  of  the  limit  state  equation  as  do  the  Taylor  series  expansions  (Rosenbluth 
1975).  They  also  require  no  assumptions  on  the  correlation  of  independent  vari¬ 
ables.  Further,  the  methods  can  be  done  on  a  calculator  or  in  a  spreadsheet  and 
do  not  require  the  laborious  calculations  of  the  Monte  Carlo  method.  Accuracy 
of  the  methods  depends  primarily  on  the  accuracy  of  the  approximate  probability 
mass  functions. 

A  common  PEM  is  a  two-point  approximation  to  the  probability  distribution, 
lumping  the  distribution  into  two  masses  at  one  standard  deviation  above  and 
below  the  mean.  This  method  is  similar  to  the  TSFD  mean  point  technique  de¬ 
scribed  above.  The  accuracy  of  the  method  is  similar  to  the  accuracy  of  the 
TSFD  method  discussed  above. 


System  Reliability 

The  reliability  of  various  modes  of  failure  or  multiple  components  in  a  design 
can  be  treated  using  a  system  analysis.  System  analyses  are  discussed  in  the 
references  and  will  only  briefly  be  discussed  herein.  System  reliability  is  deter¬ 
mined  by  combining  the  components  in  series  or  parallel  fashion,  or  a  combina¬ 
tion  of  the  two. 
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Series  system 


A  system  of  components  can  be  considered  to  be  in  series  if  failure  of  any 
one  component  leads  directly  to  system  failure.  If  an  individual  component  has 
a  probability  of  satisfactory  performance,  or  reliability,  of  Rj,  then  the  reliability 
of  a  system  of  n  components  in  series  will  be  the  product  of  the  individual 
reliabilities  as 


(26) 


Parallel  system 

A  system  can  be  considered  to  be  parallel  if  failure  of  all  components  must 
occur  for  system  failure  to  occur.  Therefore  the  reliability  of  a  parallel  system  is 
given  by 


R,  =  1  -  (1  -  R,)il  -  R^)(l  -  R^)...(l  -  /?.)...(!  -  R^)  (27) 


Parallel  series  system 

Most  projects  can  be  considered  to  be  a  collection  of  parallel  and  series  sys¬ 
tems.  The  analysis  becomes  very  complex  when  interaction  between  the  compo¬ 
nents  is  considered.  This  is  particularly  true  for  coastal  structures,  as  the  interac¬ 
tion  between  modes  of  failure  is  generally  unknown.  Several  methods  exist  for 
determining  the  reliability  of  a  complex  system.  The  most  widely  used  is  simply 
to  determine  the  upper  and  lower  bounds  using  a  separate  series  and  parallel 
analyses.  Harr  (1987)  suggests  methods  for  narrowing  the  bounds.  Another 
common  technique  for  analyzing  complex  systems  is  to  apply  weighting  factors 
to  each  component.  These  methods  are  beyond  the  scope  of  this  report  but  are 
discussed  at  length  in  the  references. 

The  most  practical  approach  may  be  to  focus  on  the  components  or  subsys¬ 
tems  that  govern  the  reliability  of  the  overall  system  (USAGE  1992).  Therefore, 
breakwater  reliability  may  be  judged  acceptable  if  the  reliability  index  for  each 
component  exceeds  a  specified  level. 
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3  Breakwater  Reliability 
Examples 


Introduction 

In  this  chapter  the  reliability  computational  methods  of  the  previous  chapter 
are  applied  to  coastal  rubble  structure  design,  and  the  methods  are  compared. 

The  reliability  indices  are  computed  for  common  limiting  states  of  performance 
including  armor  unit  stability,  armor  structural  response,  and  runup  and  overtop¬ 
ping.  Any  limit  state  equations  or  performance  functions  can  be  used  in  the 
reliability  analysis.  Those  used  herein  are  considered  accepted  practice  as  estab¬ 
lished  in  the  Shore  Protection  Manual  (SPM  1984).  In  addition  to  the  Hudson 
equation,  as  previously  described,  the  stability  equations  of  van  der  Meer  (1987) 
are  also  evaluated.  The  performance  functions  of  SPM  sections  7-II-1  and  7-II-2 
are  used  for  runup  and  overtopping,  respectively.  The  methods  of  Melby  (1989, 
1993)  are  used  for  concrete  armor  stress  evaluation. 


Hudson  Armor  Stability  Limit  State 

The  Hudson  equation  was  rewritten  in  the  form  of  a  safety  factor  (Equation 
2)  or  as  a  safety  margin  (Equations  3  and  4)  in  Chapter  2.  All  variables  in  the 
Hudson  equation  are  considered  stochastic  in  this  analysis. 

Wave  height 

The  single-storm  wave  heights  generally  follow  a  Rayleigh  probability  distri¬ 
bution  of  the  form 

p(C)=2Ce-^'  (28) 


where  C  =  H/H^,  and  is  the  root  mean  square  wave  height.  The  significant 
wave  height  is  usually  given  in  design  and  can  be  incorporated  using  the  relation 
H^=  1.416  H^.  The  mean  and  standard  deviation  required  for  the  reliability 
analysis  are  given  by  pH=0.886/frms  ^nd  o^=0.463H^^,  respectively. 
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Note  that  any  distribution  of  wave  heights  could  be  used,  including  a  mea¬ 
sured  distribution  or  an  extremal  Weibull  distribution.  Additionally,  because 
water  depth  and  wave  height  distributions  are  correlated  in  shallow  water,  any 
long-term  analysis  must  take  this  into  account.  A  number  of  techniques  exist  to 
facilitate  combination  of  the  long-term  wave  height  and  water  level  distribu¬ 
tions,  including  Monte-Carlo  simulation.  Empirical  Simulation  Technique, 
(USAGE  1996),  and  direct  integration  of  historical  distributions.  Within  this 
report,  the  wave  height  distribution  is  overly  simplified  to  focus  on  the  reliability 
calculation. 

Empirical  stability  coefficient 

The  stability  coefficient  depends  on  many  factors  including  armor  shape  and 
roughness,  wave  period,  depth,  storm  duration,  level  of  damage,  and  offshore 
bathymetry.  In  deterministic  design,  the  minimum  stability  coefficient  derived 
from  hydraulic  stability  tests  is  typically  used.  For  example,  for  two  layers  of 
rough  angular  stone  placed  on  a  structure  trunk  exposed  to  breaking  waves  ,  the 
SPM  recommends  a  minimum  stability  coefficient  of  2.  But,  as  can  be  seen  in 
Figure  5,  there  is  a  wide  variation  of  minimum  stability  coefficients,  even  for  the 
limiting  case  of  breaking  waves  on  a  breakwater  trunk.  The  variation  shown  in 
this  figure  is  only  due  to  variability  in  laboratory  mound  construction  and  vari¬ 
ability  in  wave  phasing.  Carver  (1983)  provided  data  which  suggest  a  mean  Kp 
of  2.59  and  a  standard  deviation  of  0.65  for  angular  stone  on  a  jetty  tmnk  ex¬ 
posed  to  breaking  waves.  For  concrete  armor  units,  the  variation  of  stability 
coefficients  is  equally  wide.  For  dolos,  data  of  Carver  (1983)  suggest  a  mean  of 
21.5  and  standard  deviation  of  5.6  for  similar  wave  and  structure  conditions.  For 
dolos  stability,  the  coefficient  of  variation,  defined  as  the  ratio  of  the  standard 
deviation  to  the  mean,  is  0.26,  according  to  these  data.  For  other  concrete  armor 
shapes,  a  similar  coefficient  of  variation  can  be  assumed  where  insufficient  data 
exist  to  estimate  the  value. 

Structure  seaside  slope,  armor  specific  weight,  and  armor  weight 

The  seaward  structure  slope  9  is  usually  carefully  controlled  in  construction 
and,  with  no  supporting  data,  a  coefficient  of  variation  of  0.05  is  provisionally 
used  herein.  The  specific  weight  of  stone  and  concrete  y,  is  assumed  to  have  a 
coefficient  of  variation  of  0.03  based  on  prototype  data  (Bryant  and  Mlakar 
1990).  The  final  stochastic  variable  is  the  armor  weight.  For  stone,  the  grada¬ 
tion  is  prescribed  as  75  percent  to  125  percent  and  can  be  assumed  provisionally 
to  follow  a  normal  distribution  with  a  coefficient  of  variation  of  0.10.  For  con¬ 
crete  armor,  the  coefficient  of  variation  of  weight  is  small  and  is  assumed  to  be 
0.01. 
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Iterative  TSFD  reliability  calculation 

The  reliability  calculation  using  the  ITSFD  approach  of  Equation  24  and  the 
safety  margin  Hudson  limit  state  of  Equation  3  has  been  applied  herein  using  the 
statistics  from  Table  2.  The  resulting  reliability  and  reliability  indices  versus  the 
expected  value  of  weight  of  stone  are  shown  in  Figure  6.  A  similar  plot  is 
shown  for  dolos  armor  in  Figure  7.  In  these  figures,  B  refers  to  p,  the  reliability 
index.  In  Figure  6,  the  results  for  the  equations  of  van  der  Meer  (1985)  indicated 
by  VdM  are  also  shown.  These  equations  will  be  discussed  in  the  following 
section.  As  can  be  seen  in  Figures  6  and  7,  reliability  increases  with  armor 
weight  from  zero  quickly,  approaching  R  =  1  asymptotically.  Using  the  SPM 
recommended  =  2  with  H =  1 .8  =  5.5  m  yields  a  design  armor  weight 

of  41  tons,  which  has  a  reliability  oiR~  0.98.  The  corresponding  reliability 
index  is  approximately  2.  The  reliability  index  increases  at  a  slower  rate  than  R 
with  armor  size. 


CC. 

.2 

0) 

t£. 


Hudson  B  ®  VdM  B  "  Hudson  R  VdM  R 


Figure  6.  ITSFD  solution  for  stone  stability  limit  state 
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Figure  7.  ITSFD  solution  for  dolos  stability  limit  state 


TSFD  reliability  calculation 


The  underlying  calculations  for  the  TSFD  safety  factor  approach  in  Equations 
20  and  21  with  the  partial  derivatives  evaluated  using  the  finite  difference  calcu¬ 
lation  of  Equation  23  are  shown  in  Table  3.  The  stone  weight  for  this  calculation 
corresponds  to  =  2.  It  can  be  seen  that  the  values  of  each  variable  are  all 
mean  values  except  for  a  single  perturbation  one  standard  deviation  above  the 
mean  and  one  below.  The  reliability  index  is  calculated  as  P  =  Pinp/ciinF  =  1.4. 
The  influence  factors,  given  by 


(29) 


show  the  influence  of  each  variable  on  the  overall  reliability.  It  is  clear  that,  for 
this  formulation  of  the  Hudson  performance  function,  wave  height  is  dominant 
and  that  all  other  variables  could  be  considered  to  be  constants  without  signifi¬ 
cantly  affecting  the  calculation.  This  is  due  to  the  fact  that  the  wave  height  is 
cubed  in  the  performance  function,  and  is  due  to  the  broad  wave  height  distribu¬ 
tion.  If  a  narrower  wave  height  distribution  were  used,  the  influence  of  the  sta¬ 
bility  coefficient  could  be  similar  to  that  of  the  wave  height. 
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Table  3 

TSFD  Reliability  Calculation  for  Hudson  Limit  State 

Armor 

Weight 

W 

kg 

Struc¬ 
ture 
Slope 
cot  6 

Hudson 

Coeff. 

KD 

Armor 

Specific 

Weight 

Y 

kg/m' 

Wave 

Height 

H 

m 

Safety 

Factor 

F 

Finite 

Differ¬ 

ence 

dF/dX 

infiu- 

ence 

Factor 

a^/Ia' 

% 

45013 

2 

2.59 

2,405 

2.70 

10.97 

0.1003 

0.1 

36829 

2 

2.59 

2,405 

2.70 

10.75 

40921 

2.1 

2.59 

2,405 

2.70 

11.40 

0.0500 

0.0 

40921 

1.9 

2.59 

2,405 

2.70 

10.32 

40921 

2 

3.24 

2,405 

2.70 

13.58 

0.2564 

0.3 

40921 

2 

1.94 

2,405 

2.70 

8.13 

40921 

2 

2.59 

2,477 

2.70 

12.18 

0.1271 

0.1 

40921 

2 

2.59 

2,333 

2.70 

9.62 

40921 

2 

2.59 

2,405 

4.11 

3.08 

-1.740 

99.5 

40921 

2 

2.59 

2,405 

1.25 

99.78 

Comparison  of  reliability  methods  for  Hudson  performance  function 

Results  from  the  Level  II  calculations,  applied  to  the  Hudson  limit  state,  are 
summarized  in  Table  4  for  a  single  stone  eumor  weight  and  other  parameters 
listed  in  Table  2.  The  stone  weight  of  41  tons  utilized  corresponds  to  Kd=  2,  as 
the  most  conservative  recommendation  in  the  SPM.  It  is  clear  from  Table  4  that 
the  TSFD  gives  a  more  conservative  result  than  the  ITSFD  method.  The  TSFD 
is  simpler  but  accuracy  is  sacrificed.  Note  that  reliability  was  computed  using 
the  ITSFD  approach  with  two  different  versions  of  the  Hudson  equation  perfor¬ 
mance  function.  The  resulting  reliability  and  reliability  index  were  identical. 
Therefore,  the  ITSFD  was  verified  to  be  invariant  to  the  form  of  the  limit  state. 


Table  4 

Comparison  of  Taylor  Series  Approximations  of  Reliability  for 
Hudson  Limit  State 

Reliabiiity  Computation 
Method 

Reliability 

R 

Reliability 

Index 

P 

Difference  in  R 
From  ITSFD 

% 

ITSFD  (Equation  3) 

0.98 

2.0 

- 

ITSFD  (Equation  4) 

0.98 

2.0 

- 

TSFD 

0.91 

1.4 

-7 
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Van  der  Meer  Armor  Stability  Limit  State 


Armor  stability  equations  for  two-dimensional  trunk  sections  were  proposed 
by  van  der  Meer  (1987)  for  stone,  cubes,  Accropodes,  and  Tetrapods.  Although 
not  verified  for  general  use  in  three-dimensional  breakwater  applications,  these 
equations  provide  insight  into  various  parameters  of  interest  in  stability,  includ¬ 
ing  wave  period,  porosity,  degree  of  damage,  and  duration  of  storm,  that  are  not 
included  explicitly  in  the  Hudson  equation.  The  equations  have  been  rearranged 
into  limit  state  relations  as  follows: 

STONE,  ?  < 


g  =  K^i6.2AD„  (30) 

STONE,  ?  > 

g  =  (31) 


where  ^  =  ttin0/(5„)  ^,  =  H^/L„,  L  =  gT^Hn,  is  the  mean  wave  period,  P  is 

the  structure  porosity,  S  is  the  degree  of  damage,  N  is  the  number  of  waves,  and 
K^,  and  are  empirical  coefficients. 


CUBES 


5  S. 


-0.7 


‘  04  ^ 

6.7—  +  1.0 


iV 


0.3 


-  H. 


) 


(33) 


TETRAPODS 


S  =  s, 


-0.2 


3.75—2—  +  0.85 
^0^ 


-  H. 


(34) 
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where  is  the  zero-damage  number  of  waves.  Table  5  provides  sample  statis¬ 
tics  for  application  of  the  various  parameters  in  the  equations  listed  above. 

Figure  6  compares  the  ITSFD  method  applied  to  the  van  der  Meer  equations 
to  that  from  the  Hudson  equation.  For  these  input  conditions,  the  two  equations 
yield  very  similar  reliability  values  for  the  range  of  armor  weights.  Comparing 
van  der  Meer  equations  with  the  Hudson  results,  it  is  clear  that  the  van  der  Meer 
equations  yield  a  higher  reliability  than  does  the  Hudson  equation. 


Table  5 

Statistics  for  Van  der  Meer  Equation  Limit  State  Variabies 

Random 

Variable 

Mean 

Standard 

Deviation 

Coefficient  of  Varia¬ 
tion 

M 

o 

c 

Stone  nom. 
diameter,  D„ 

Range 

Range 

0.01 

Stone  stab.  eq.  coef¬ 
ficient,  Ks, 

6.2 

0.4 

0.065 

Stone  stab.  eq.  coef¬ 
ficient. 

1.0 

0.08 

0.08 

Stone  specific  grav¬ 
ity,  Sr 

2.6 

0.08 

0.03 

Tetrapod  nom.  diam¬ 
eter,  D„ 

Range 

Range 

0.0 

Tetrapod  eq.  coeffi¬ 
cient,  K, 

1.0 

0.1 

0.10 

Concrete  spec,  grav¬ 
ity,  Sr 

2.34 

0.07 

0.03 

Cube  nom. 
diameter, 

Range 

Range 

0.0 

Cube  eq.  coefficient, 

Kc 

1.0 

0.1 

0.10 

Wave  height,  H 

2.70  m 

1.41  m 

0.52 

Wave  steepness,  S„ 

0.04 

0.01 

0.25 

Porosity,  P 

0.5 

0.05 

0.10 

Number  of  waves,  N 

2000 

600 

0.03 

Structure  slope,  cot9 

2.0 

0.1 

0.05 

Concrete  Armor  Structural  Limit  State 

The  concrete  armor  unit  strength  design  methods  of  Melby  (1989, 1993) 
provide  a  Level  HI  reliability-based  analysis  method  where  various  loading  dis- 
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tributions  are  convolved  to  determine  a  design  load  distribution.  Here,  the  de¬ 
mand  is  characterized  as  a  design  tensile  stress  for  the  armor  layer  and  the  ca¬ 
pacity  is  the  tensile  strength  of  the  unreinforced  concrete.  A  design  probability 
of  exceedance  is  specified  in  order  to  determine  the  design  stress  from  the  total 
stress  probability  distribution.  This  design  stress  is  compared  with  a  predefined 
capacity,  a  fatigue-reduced  concrete  strength.  For  this  Level  III  analysis 
method,  the  design  probability  of  exceedance  E  is  used  to  determine  a  design 
stress  from  the  maximum  stress  design  probability  distribution.  This  design 
exceedance  is  characterized  as  the  probability  of  exceeding  the  design  stress  in 
the  armor  layer. 

To  relate  the  Level  III  methods  to  the  Level  II  analysis  method,  the  design 
exceedance  is  related  to  P  through  the  basic  reliability  relation  E  =1-R  =  1-$(P). 
For  concrete  armor  units,  the  limiting  state  of  structural  performance  is  given  by 

=  0  (35) 


where/,  is  the  flexural  tensile  strength  and  is  the  maximum  principal  tensile 
stress  in  the  armor  unit.  Ellingwood  et  al.  (1980)  summarized  a  number  of  stud¬ 
ies  of  the  distribution  of  the  tensile  strength  of  concrete  and,  based  on  this  work, 
a  coefficient  of  variation  of  18  percent  is  suggested  (Mlakar  1995).  The  distri¬ 
bution  of  Oj  is  computed  in  PC-ARMOR  (Melby  1989, 1993)  and  used  in  the 
reliability  analysis  to  determine  the  reliability  with  respect  to  armor  strength. 

An  example  calculation  was  done  using  PCARMOR  to  compute  a  mean  of 
3.05  MPa  and  a  standard  deviation  of  1.46  MPa.  The  results  of  a  simple  analysis 
are  shown  in  Figure  8. 


tx: 


o 

a: 


Armor  Strength  (MPa) 


strength  R 


Strength  B 


Figure  8.  ITSFD  solution  for  dolos  structural  limit  state 
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Runup  and  Overtopping  Limit  States 


As  stated  previously,  the  performance  functions  of  SPM  sections  7-II-1  and 
7-II-2  are  used  for  runup  and  overtopping,  respectively.  The  equation  for  runup 
is  given  as 


=  1[1  -exp(-^2)] 

Urn 


(36) 


where  R  is  the  significant  runup,  H,  the  significant  wave  height,  and  ^  the  surf 
similarity  parameter  as  used  previously.  A  safety  margin  limit  state  can  be  de¬ 
fined  as  the  difference  between  the  freeboard/and  the  runup 


g=f-R=f  -  l.m[l  -exp(-^2)] 


The  safety  factor  limit  state  could  be  stated  as  F=f/R. 
For  overtopping,  the  equation  given  in  the  SPM  is 

Kjjj  =  (0.51  -  O.llB/AXl.O  -f/R) 


(37) 


(38) 


where  Kj-o  is  the  coefficient  for  transmission  by  overtopping,  B  the  crest  width, 
and  h  the  crest  height  from  the  bottom.  The  safety  margin  limit  state  equation 
can  be  expressed  as 


8 


allow 


K^H 


(39) 


Both  runup  and  overtopping  have  been  incorporated  into  the  ITSFD  reliabil¬ 
ity  methods  and  the  safety  margin  relations  encoded  similar  to  the  example 
given  in  Appendix  A.  Input  statistics  for  the  calculations  are  shown  in  Table  6. 
Results  from  an  example  computation  are  shown  in  Figures  9  and  10. 
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Reliability  Index, 


0.6  1.2  1.8  2.43.1  3.7  4.3  4.9  5.5  6.1 
Freeboard  (m) 


Overtop  R  *  Overtop  B 


Figure  10.  ITSFD  solution  for  overtopping  limit  state 
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Reliability,  R 


4  Conclusions 


In  this  report,  several  level  II  reliability  techniques  are  discussed  and  applied 
to  breakwater  design.  Using  these  methods,  the  reliability  and  reliability  index 
are  determined  for  the  dominant  performance  functions  of  a  breakwater  includ¬ 
ing  stone  and  concrete  armor  stability,  concrete  armor  structural  response,  runup, 
and  overtopping.  The  methods  include  Taylor  series  finite  difference  (TSFD) 
methods,  which  are  shown  to  be  easily  computed  using  a  handheld  calculator  or 
spreadsheet.  The  TSFD  method  is  shown  to  yield  reasonable  accuracy  for  pre¬ 
liminary  comparison  of  various  alternatives.  An  invariant  iterative  TSFD 
method  and  associated  FORTRAN  program  are  shown,  which  provide  an  im¬ 
proved  approximation  to  the  reliability  for  design  and  evaluation  purposes. 
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Appendix  A  Computer  Pro¬ 
gram  Listings 


Appendix  A  Computer  Program  Listings 


A1 


PROGRAM  RELIABLE 


Level  II  iterative  invariant  first  order  Taylor  series  design  point 
reliability  engine.  Program  computes  reliability  and  reliability  index 
for  a  range  of  independent  variable  values.  The  algorithm  is  invariant 
with  respect  to  the  form  of  the  limit  state  equation. 


c 

The  limit  state  equation  is  computed  in  the  subroutine  LIMSTATE. 

c 

The  program  computes  various  values  of  the  independent  variable  in 

c 

o 

subroutine  DEP_MEAN. 

u 

c 

Data  are  input  from  a  user-edited  file  called  R  INPUT.DAT  and  data  are 

c 

Output  to  a  file  called  R_OUTPUT.DAT. 

u 

c 

Melby  11/95 

c 

SCALAR  DEFINITION 

c 

betaz 

Reliability  Index  at  z 

c 

betaz_p 

Reliability  Index  at  design  point 

c 

g 

Limit  state 

c 

gz 

Limit  state  at  z 

c 

I.J 

Indices  for  random  variable 

c 

ik# 

Iteration  kounter  # 

c 

iki 

Index 

c 

LAMBDA 

Linear  coefficient  between  z’  and  g 

c 

NumRV 

Number  of  random  variables 

c 

R 

Reliability 

c 

cov 

std  /  mean 

c 

Q 

zd 

z  at  design  point 

c 

VECTOR  DEFINITION 

c 

grad_g{NumRV) 

Gradient  of  limit  state 

c 

Mean(NumRV),m 

Means 

c 

Std(NumRV),s 

Standard  devs 

c 

X(NumRV) 

Random  variables 

c 

z(NumRV) 

Standard  normal  random  variables 

c 

zp(NumRV) 

Design  point 

************************* 


PARAMETER  (NRV=10) 

REAL  mean(NRV),  Std(NRV),  z(NRV),  zp(NRV).  grad_g(NRV) 
REAL  Const(NRV),  cov{NR\/),  x(NRV),  lambda 
CHARACTER  textc(NRV)*15,  textR(NRV)*15,  header*60,  texth*15 


Open  I/O  files 


OPEN(5,  FILE='RJnput.dat',  status='old') 
OPEN(2,  FILE=’R_output.dat’,  status=' unknown’) 


Read  input  data  file 


PRINT*,’  Program  will  now  read  input  file  named  R_input.dat’ 
PRINT* 


PRINT*,’  Input  statistics  for  each  variable  are  input  in  this  file.  For  each  random’ 
PRINT*,’  variable,  the  mean,  std.  Dev.,  and  coef.  of  variation  are  read.  Only  the  mean’ 
PRINT*,’  and  s.d.  are  required  for  each  variable  except  the  last  one  input,  which’ 
PRINT*,’  will  be  treated  as  the  dependent  variable  and  the  reliability  computed  over  a’ 
PRINT*,’  range  of  values.’ 

PRINT* 

PAUSE 

READ(1,’{A60)’)header 

PRINT*,header 

READ(1,’(A60)’)header 

PRINT*,header 

READ(1  ,*)texth,  numRV,  NumConst,  NumDepVar 
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PRINT* 

PRINT*,’ - Values  read  from  input  file - ’ 

PRINT*,’  Number  of  random  variables’, NumRV 
PRINT*,’  Number  of  constants’, NumConst 
PRINT*,’  Number  of  variable  calculatlons’,NumDepVar 
PRINT* 

READ(1,’(A60)’)header 


PRINT*,’ - Random  variable  statistics  read  from  input  file 

PRINT*,’ Param  Mean  Std  Dev  CoefofVar’ 
PRINT*,’ - ’ 


DO  5,  ikO  =  1 ,  NumRV 

READ(1  ,*)  textrv(lkO),  mean(ikO),  Std(ikO),  Cov(ikO) 

PRINT*,  textrv(ikO),  mean(ikO),  Std(ikO),  Cov(ikO) 

5  CONTINUE 
READ(1,’(A60)’) 

PRINT* 

PRINT*,’ - Constants  read  from  input  file - ’ 

PRINT*,’  Constant  Value  ’ 

PRINT*,’ - ’ 

DO  6,  ikO  =  1 ,  NumConst 
READ(1  ,*)  textc(ikO),  Const(ikO) 

PRINT*,  textc(ikO),  Const(ikO) 

6  CONTINUE 
PRINT* 

PRINT*,’  Values  OK?  If  not,  press  CTRL-C  and  modify  the  input  file’ 
PAUSE 


c  Write  output  file  headers 


WRITE(2,*)’  Output  data  file  from  RELIABLE’ 

WRITE(2,*) 

WRITE(2,*)’ - Constants  read  from  input  file - ’ 

WRITE(2,*)’  Number  of  random  variables,  constants,  and  dependent  var.  calcs.’ 

WRITE(2,*)’ - ’ 

WRITE(2,*)NumRV,  NumConst,  NumDepVar 
WRITE(2,*) 

WRITE(2,*)’ - Parameters  read  from  input  file - ’ 

WRITE(2,*)’  texth  =  ‘,texlh,’  NumRV  =  ’,  numRV,  ‘  NumConst  =  ’,  NumConst 
WRITE(2,*) 

WRITE{2,*)’  Parameter  Mean  Std.  Dev.  Coef.  of  Var.’ 

WRITE(2,*)’ - ’ 

D0  7ikO  =  1,  NumRV 

WRITE(2,*)  textrv(ikO),  mean{ikO),  Std(ikO),  Cov(ikO) 

7  CONTINUE 
WRITE(2,*) 

WRITE(2,*)’  Constant  Value’ 

WRITE(2,*)’ - ’ 

DO  8  ik0=1 ,  NumConst 
WRITE(2,’)  textc(ikO),  Const(ikO) 

8  CONTINUE 
WRITE{2,*) 

WRITE(2,*)’ - Output  Results - ’ 

WRITE(2,*) 

WRITE(2,*)’  Mean  Beta  R’ 

WRITE{2,*)’ - ’ 

c 

c - 

c  Initialize  variables 


DO  to  kount  =  1  ,NumRV-1 
z(kount)  =  0.0 
10  CONTINUE 
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A3 


c - 

c  Begin  Computations 

c - 

DO  20  ikl=1 ,  NumDepVar 

CALL  DEP_MEAN(ikl,  mean,  NRV,  Const,  NumRV) 

Std(NumRV)  =  cov(NumRV)  *  Mean(NumRV) 

Z(NumRV)  =  0.0 
betaz  =  0.0 
30  CONTINUE 

CALL  LIMSTATE(z,  x,  std,  mean.  Const,  g,  NRV,  NumRV) 

gz  =  g 

grad_g0  =  0.0 
zd  =  0.0 
c 

DO40ik2  =  1,NumRV 
z(ik2)  =  z(ik2)  +  0.001 

CALL  LIMSTATE(z,  x,  std,  mean.  Const,  g,  NRV,  NumRV) 

z(ik2)  =  z{ik2)  -  0.001 

grad_g(ik2)  =  (g  -  gz)/0.001 

grad_g0  =  grad_g0  +  grad_g(ik2)’grad_g(ik2) 

zd  =  zd  +  z(ik2)  *  grad_g(ik2) 

40  CONTINUE 
c 

lambda  =  (zd  -  gz)  /  grad_g0 

0 

c  Find  zp  and  betaz_P 

betaz_P  =  0.0 
DO  50  ik3  =  1, NumRV 
zp{ik3)  =  lambda  *  grad_g(ik3) 
betaz_P  =  betaz_P  +  zp(ik3)  *  zp(ik3) 

50  CONTINUE 

betaz.P  =  SQRT(betaz_P) 

IF(ABS(betaz_P  -  betaz)  .LT.  0.001)  GOTO  90 
betaz  =  betaz_P 
c 

DO  60  ik4  =  1, NumRV 
z(ik4)  =  zp(ik4) 

60  CONTINUE 

GOTO  30 

90  CONTINUE  !  solution  converged 

c 

c  Hastings  rational  approximation  for  normal  reliability 

R  =  1.  -  0.5/(1.  +  0.196854  ’betaz_P  +0.115194*betaz_P*betaz_P+ 
&  0.000344*betaz_P**3  +  0.01 9527*betaz_P**4)**4 

c 

WRITE{2,100)  mean(NumRV),  betaz_P,  R 
100  FORMAT(1x,  3(F15.3)) 

20  CONTINUE 
CLOSE(I) 

CLOSE(2) 

STOP 

END 
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SUBROUTINE  limstate(z,  x,  s,  m,  C,  g,  NRV,  NumRV) 
REAL  z(NRV),  x(NRV),  s(NRV),  m(NRV),  C(NRV) 
c  Hudson  Stability  Equation 
c  Random  Variabies,  mean,  standard  deviation 


c  Conor  Specific  Wt 

Gamma_r 

m(1) 

s(1) 

c  Structure  Slope 

cot  theta 

m(2) 

s(2) 

c  Coefficient 

KD 

m(3) 

s(3) 

c  Wave  Height 

H 

m(4) 

s(4) 

c  Armor  Weight 

W 

m(5) 

s(5) 

c  Constants 

c  H20  Specific  Weight 

Gamma_w 

c(1) 

c 

DO  10j=1, NumRV 
xO)  =  z{j)  *  sG)  +  mC) 

10  CONTINUE 

Delta  =  x(1)/C(1)-1. 

g  =  x(2)  *  x(3)  *  Delta**3  *  x(5)  -  x(1)  *  x(4)**3 

RETURN 

END 


SUBROUTINE  dep_mean(kount,  m,  NRV,  C,  NumRV) 
REAL  m(NRV),  C(NRV) 
c  Hudson  Stability  Equation 
c  Random  Variables,  mean,  standard  deviation 


c  Conor  Specific  Wt 

Gamma_r 

m(1) 

s(1) 

c  Structure  Slope 

cot  theta 

m(2) 

s(2) 

c  Coefficient 

KD 

m{3) 

s(3) 

c  Wave  Height 

H 

m(4) 

s(4) 

c  Armor  Weight 

W 

m(5) 

s(5) 

V 

c  Constants 

c  H20  Specific  Weight 

Gamma_w 

c{1) 

c 

Delta  =  m(1)/C(1)-1 

m(NumRV)  =  FLOAT(kount)*  m(1)*  (m(4)/Delta)**3  /  m(3)/  m(2) 

RETURN 

END 
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Hudson  Damage  Function  Reliability  Input  in  File  RJNPUT.DAT 
Sample  Stone  Stability  Run,  Melby,  19  Nov  1995 
Constants  5  1  11 

Input  Random  Variable  Statistics 


Gamma_r 

2405.0 

72.0 

-1. 

cot 

2. 

0.1 

-1. 

KD 

2.59 

0.65 

-1. 

H 

2.70 

1.41  ■ 

-1. 

W 

-1.0 

-1.0 

0.2 

Input  Constants  in 

Limit  State 

Equation 

Gamma_w 

1024.0 

Output  data  file  from  RELIABLE 

- Constants  read  from  input  file - 

Number  of  random  variables,  constants,  and  dependent  var.  calcs. 


5  1  11 

- Params  read  from  input  file - 

texth  =  Constants  NumRV  =  5  NumConst  =  1 

Parameter  Mean  Std  Dev  Coef  of  Var 


Gamma_r 

2405.0 

72.000 

-1.0000 

cot 

2.0000 

0.1000 

-1.0000 

KD 

2.5900 

0.6500 

-1.0000 

H 

2.7000 

1.4100 

-1.0000 

W 

-1.000 

-1.000 

0.20000 

Constant 

Value 

Gamma_w 

1024.0 

Output  Results 


Mean 

Beta 

R 

3768.453 

0.000 

0.500 

7536.913 

0.483 

0.686 

11305.357 

0.815 

0.792 

15073.811 

1.074 

0.858 

18842.263 

1.289 

0.901 

22610.715 

1.473 

0.930 

26379.168 

1.635 

0.949 

30147.622 

1.780 

0.963 

33916.074 

1.912 

0.972 

37684.526 

2.032 

0.979 

41452.977 

2.142 

0.984 
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